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We investigate the change in the decay rate of a muon caused by embedding it in the field of a 
laser. A previous paper found that the change could be large, as much as an order of magnitude. 
] We find the more intuitive result that the change is small and give analytic expressions for the small 
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corrections. 

PACS numbers: 13.35.Bv, 13.40.Ks, 14.60.Ef, 42.62.-b 
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There is interest and work on the properties of elementary systems when they are placed in strong electromag- 
netic fields [lj, [2], |3j. Recently, in an attempt to extend this work to unstable systems, Liu, Li, and Berakdarjj] 
(LLB) calculated the effect of a strong laser field on the decay rate of muons. They found the lifetime could be 
changed from its normal value of 2.2 x 10 6 seconds to 5 x 10 7 seconds or even less. This conclusion was chal- 
lenged by Narozhny and Fedotov, who offer an abbreviated calculation to support their criticism [1, 0| . If LLB 
were correct, this would be a very interesting result. We have done our own calculation and, unfortunately, also 
reach a very different conclusion. We find the effects of the laser to be very small and give explicit expressions 
for these small effects. 

£■ — ' Although the idea of the problem is straightforward, the actual calculation is somewhat complicated and LLB 

^sO . did the complicated part numerically. We do everything analytically. Because of this difference, and because we 
get such a different result, we will present our calculation in some detail and only after our results are apparent 
will we compare with LLB. Also because our calculation is analytic, we don't need to make definite choices 
Q^ ■ about the properties of the laser. We will assume only that the energies of the photons are about 0.1 — 1.0 eV 
and that the magnitude of the laser field amplitude is 10 6 — 10 7 V/cm (as used in LLB). The next section gives 
our work, the following section compares with, and discusses, LLB, including the fact that we have somewhat 
different starting points; we use the Volkov wavefunction 0, HI, U flCj | for circular polarization while they use an 
approximation to the wavefunction for linear polarization. The last section repeats our conclusions. 

X 

H 

2. FORMALISM 

The process is muon decay, 

H~(P) — ► e~(p) + u e {qi) + i/ M (<fc) , (1) 

where the arguments are our labels for the momentum. We will assume the photons from the laser are along the 
z-axis with momentum = (u>, 0, 0, ui) and circular polarization, 

A p '{x) = arii cos k ' x + an 2 sni k ■ x , (2) 
< = (0,1,0,0), (3) 
n» = (0,0,1,0). (4) 
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The electron wave function is then 

Mx) = e"' 1 ^"" sin k - x e~ l ^ k - x ~ lp - x ( 1 + -^L) u (p) (5) 

V 2p ■ k / 

where we have taken the electron to be in the xz plane and thereby avoided a factor —i^^Py cos k ■ x in the 
exponential. From the second exponential factor the electron has an effective momentum and mass 

t = ^ + (6) 

2 2,2 2 /■r 7 \ 

m = m + e a , (7) 

where mo = 0.511 MeV. Note that q ■ k = p ■ k and q x = p x - Following the standard proceedure we use the 
generating function for Bessel functions [111]. 



E t'Mz) (8) 



to rewrite the first factor in ([5]) as 



E Je(D)e- uk - x (9) 



with Z? = S££k Momentum conservation is then 

p- k 

Pf+tk" = q^ + qf + q^ (10) 
and the matrix element, for a given value of £, is 

Re = -^S(g 2 )7 Q (l-75)w(P)u(p)[Ao + Ai^i^ + A 2 ^ 2 ^7 Q (l-7 5 )?;(gi) (11) 

where 

A = J e (D) (12) 

Ai = ^(%(D) + Jn(5)) (13) 

A 2 = y^(J m (i?)-J.-i(i))). (14) 

Note that the argument of the electron spinor is still p. 

We square the matrix element in the usual way, using FORM|ll|, and integrate out the momentum of the 
neutrinos in the usual way using 

/ ^i^i si{Q - qi - q2)q ^ = ^(Q 2 -9 Q/3 + 2g^)e(Q 2 ). (is) 

The total width is then 

oo 

r = £ r ^ ( 16 ) 

£=-oo 

with the width for each £ given as an integral over the electron energy and angle 
where E is q° (the Jacobian connecting p^ and q^ is unity) and M is the muon mass. 
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The square of the matrix element, summed over spin and integrated over the neutrino momentum, is 
|T £ | 2 = 128G 2 {j 2 (L>)[3P- g(M 2 + to 2 ) - 4(P • q) 2 - 2M 2 m 2 ] 

2 2 

+ jj{D)^—[2q ■ k(M 2 -p.g)+p. fc(4P • q - 3M 2 - to 2 )] 
zq • /c 

+ £ J 2 {D) [q ■ k{2M 2 - 4P • q) - P ■ k(3M 2 - 8P ■ q + 3m 2 )] 
e 2 a 2 

+ 2£J 2 AD) -P ■ k(q- k- P ■ k) 

q ■ k 



A£ 2 Jj{D)P ■ k(q- k- P ■ k) 

-\J 2 i+i{D) + Je-t{D)] [ - 4P • qP- k + 3M 2 P • k + m 2 P ■ k + 2q- k(P ■ q- M 2 )] 



e 2 a 2 



4q - k l 

2 2 

+ ^nji+m + jump ■ k(p-k~ q -k) 

+ ie(Pq,k,n 2 )-^-rJ e (D)[J e+1 (D)- J^D^M 2 + m 2 - AP ■ q) 
Zq- k 

pa 

+ ie(P,q,k,n 2 ) -I J t (D)[J i+1 (D) - J / _ 1 (D)](2P • k-q-k) 

q ■ k 

2 2 

+ ie{P,k,n x ,n 2 ) -^—{J 2 l+1 {D) - jj^D^M 2 + m 2 - AP ■ q) 
Aq ■ k 

2 2 

+ ie(P,fc,»i s n 2 )^r^[^ +1 (D)- Jf_ 1 (D)](2P- k ~2q ■ k) 
2q ■ k 

2 2 

+ ie{k,q,n x ,n 2 ) ^-[jf +1 (D) - j}_ x {D)]{M 2 - P ■ q)} (18) 
where e 0123 = — {. 

If we set a — then £> = so J £ = for £ ^ and Jo(0) = 1 and we get the usual expression for muon decay 
from the first line of (JT5J) 

r° = ^Ml (19) 

192tt 3 1 J 

where terms proportional to the electron mass have been neglected [15]. Note that r° is not the same as r £= o- 

The crucial thing is the limits on the integrals in Eq. 07)) . These are determined by the 8 function, 

, ,f ,i f M.+eu ,i 

dE / dz<d(Q 2 ) = I dE I dz + dE dz, (£ > 0) (20) 



<!E I d-A-)(Q 2 ) = I dE J dz + J dE J dz, (£ < 0) (21) 



where 



_ M 2 + 2M^-2£(M + M 
ZL(S) - 2£jE ' (22) 

It is important to notice that Zl(^-) = —1 and zl(*4- + £u>) = +1. We can rewrite these limits, for both signs 
of I, as 



J dE J dzQ(Q 2 ) = Jj dE J dz + J J dE J ^ 
We will treat these two terms separately. 



dz . (23) 

zl(E) 



Consider the first term in (|23[) . Because the limits of integration are not functions of t we can immediately do 
the sum over I in Eq. (fT6|) . The only things we need to do the sums are the relations 



E Je(z)Je+n(z) = J n (0), (24) 

t=-oo 

which follows from 9.1.75 in Abramowitz and Stegun(l4|, the recurrsion relation for Bessel functions, £ Je(z) = 
§(Je+i(z) + J(-i(z)), and J-g(z) = (-iyjg(z). The most important sum is 

oo 

E J e( D ) = 1 ( 25 ) 

e=-oo 

because that replaces the J 2 (D) in the first line of (fT5|) by unity. After doing all the sums Eq. (fT5|) becomes 

oo 

E \ T 4 2 = 128 G 2 {3P ■ q(M 2 + to 2 ) - 4(P • q) 2 - 2M 2 m 2 

l=-oo 

e(P,k,n u n 2 ){2P- k - 2q ■ k)} (26) 



.e 2 a 2 



q ■ k 

where there has been a lot of cancellation. This can be easily integrated by hand 

r= Wi + 8 ^-inM_^ ln M + ^n (27) 

V M 2 1 3 to M to 2M)J y ' 

where again the first line of ([2^]) gives T . Since ea is much less than 1 MeV the additional laser dependent terms 
are very small, smaller than the electron mass terms which were neglected in (|19[) . The terms which are even 
further suppressed by the factor of ui/M come from the e term in (|26|) . 

Now consider the correction from the second term in (|23[) . For ea ~ 10~ 4 MeV and uj ~ 1 eV, the argument of 
the Bessel function, Z?, varies with E and z from zero to ~ 1.5 x 10 4 . Bessel functions become very small once 
the index, £, becomes greater than the argument so £ is limited by £ < D . In other words lo£ will always be 
much less than 1 MeV and this correction will be small because the range of the energy integration is small. So 
define a function of the energy as 

F{E,lu,£) = J dEE(3 J dz\T}\ 2 (28) 
where \Tg\ 2 is given by Eq. (jTSJ) without the prefactor of 128 G 2 . Then the correction to the width is given by 

T^-^E [F{E = -+l^u,t)-F{E=-^,£)]. (29) 

f=-oo 

Now if we expand the first term in a Taylor series 

M M s dF , £ 2 u 2 d 2 F, , . 

F(E= — +tu>,u>,£) = F(E=—,u>,t)+e U —\ E= q + —-^\ E= q.+--- (30) 

then, since z L {^) = —1, we can again do the sum over £ before we do the integral. We will approximate (f2"9"| 
by keeping only the first nonzero term; since 
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the correction, Eq. 



is 



r° 



4w 
Af5 



dz V £\T e \ 2 



(32) 



Now in Eq. (|18|) the only nonzero sums come from lines 3, 4, 7, 9, 10, and 12 and after a bit of work we get 

M 



r° 



8- 



M 2 



{( 



1 



M 



M 2 



In- 



1 



M 



M 2 J 



(33) 



These terms are small but not necessarily smaller than those in (f2"7| . In fact the largest term cancels and the 
total effect of the laser, the sum of (|27l) and (f3"3")> . is to change the width from T° to 



i 



r u n +8 



M 2 



(M +4 ^ ln m 



3 w 
2M 



>/ 2 



(34) 



Let's be clear about what has been done. Eq. (|2"7| is an exact solution for the part of the phase space where 
the integral over the electron energy is between m and M- . Eq. (|33p is an approximation for the part of the phase 
space where the integral is between M- and ^-+£u>, in that it is the first nonzero term in an expansion in powers 
of £oj. As a consistancy check we have calculated the contribution to the second nonzero term in the expansion 
(the second derivative term in (|3H)) ) from the first line of (fTB|) (which should be the largest contribution) and 
found 



r c „ e 2 a 2 m 



M 

- -1 + ln — 

m 



(35) 



which is the same as zero since we did not keep terms ~ jj in (|27p or (f3"3"|) 



3. COMPARISON WITH LLB 



We have assumed the laser radiation is circularly polarized because, as is well known[2|, that case is easier. 
For linear polarization the wavefunction in Eq. ([5]) is multiplied by an extra factor 

and the second term in the definition of Eq. ©, has an extra \. Because of this term the generating function 
for Bessel functions must be used twice and the square of the matrix element involves two sums over Bessel 
functions. LLB use linear polarization but approximate the wavefunction by omitting this term thus avoiding 
the double sum. Dropping this term is not the same as doing circular polarization because the vector potential, 
is different. The effect of this approximation on the muon lifetime is unknown but it seems unlikely this 
could explain the huge difference in our results. Given their definition of the wavefunction we agree with their 
expression for the S matrix. We do not, however, agree with their expression for the partial width. We believe 
their We = / dE ■ ■ ■ should be We = 48 ff M / dEE (We is what we call IV) Otherwise their L° would 
be too big by a factor of |. This is probably a misprint. If not, it would partially explain their small lifetimes. 

LLB find a shorter lifetime which means they find a larger width. On the other hand they don't get much 
contribution from what we call negative I. (We have defined £ differently - what we call negative £ LLB call 
positive.) If we consider only positive £ then, for example, Eq. (|2"5|) is replaced by 

= -(l + J 2 Q {z)) (37) 

1=0 

and r°, the contribution from the first line of (| 18[) . is replaced by approximately ^L°. But this is a change in 
the wrong direction. LLB get a larger width despite having little contribution from negative I. As we have seen 
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the other terms in (fT8"|) go as e 2 a 2 /M 2 . Furthermore a longer lifetime might be easier to understand as a kind 
of Zeno effect; the interactions with the laser photons make the decay keep starting over. But that is not what 
LLB find. 

We disagree with LLB on the lower limit of the electron energy. We get the lower limit to be the (effective) 
electron mass, m, they get m + iu. Again it is hard to believe that makes much difference. 

So about the only place left to look for a difference is the numerical integrations. To do these requires definite 
values for ea and u>. We used ea — 1.69 x 10~ 4 MeV and u> = 1.17eV (the Nd:YAG laser of LLB) and integrated 
(fT8|) using the limits (|20|) or (j2"Tj) . A few results are shown in the table 



£ 


ftr 





3.463 


10" 


-3 


1 


3.455 


io- 


-3 


10 


3.427 


io- 


-3 


100 


1.921 


io- 


-3 


200 


6.940 


io- 


-4 


500 


7.391 


io- 


-5 


1000 


1.013 


io- 


-5 


5000 


7.866 


10" 


-8 



We used routines for the Bessel functions from Numerical Recipes fl2l] and the integration routine VEGAS[13]. 
The partial widths for negative £ were indistingushable from those of positive £ for a given \£\. The integration 
was sufficiently fast that we could do each of the for £ up to 500 with the result 

500 

^ J2 = 0.96 (38) 

£=-500 

where we estimate the error from the numerical integration to be less than 0.01. If we assume the £ dependence 
of is linear from £ = 500 to 1000 we get another contribution to (|38|) of 0.04; if we make the same assumption 
for I from 1000 to 5000 we get another 0.04. These are surely overestimates because the £ dependence must fall 
faster than linear. We could do a better job for \£\ > 500 but these are sufficient to show the magnitude of the 
contribution to the total width that could be expected from higher £. Thus the conclusion is that the total width 
cannot be very different than r°, in agreement with our more precise arguments above. 



4. CONCLUSIONS 

We have considered muon decay in the electromagnetic field of a laser. Our discussion is entirely analytic with 
the only approximation a Taylor series expansion of the squared matrix element for the region of electron energy 
between 4f and 4f + £uj, a distance of less than 10~ 2 MeV. We find the effect of a laser on muon decay is very 
small, of order e 2 a 2 /M 2 , e 2 a 2 w/M 3 , or e 2 a 2 u 2 /M 4 , where ea ~ 2 x 10~ 4 MeV, uj ~ leV and M is the muon 
mass, 105.66 MeV. Given the dimensions of the Fermi coupling constant the decay width must have five powers 
of energy. Once we find the coefficient of the M 5 term is the same as in the absence of the laser we know the 
effect of the laser is very slight because the only other energies in the problem are ea, u>, and m. We dropped the 
electron mass where possible so there are surely corrections of order e 2 a 2 m 2 /M which are numerically larger 
than some of the terms we included. But all of these corrections are tiny, smaller than the known corrections of 
order m 2 /M 2 . Our result is given by Eq. (f34|) . 
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